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^ ' Abstract. We review the significance determination for emission lines in the Chandra HETGS 

spectrum for GRB 0208 13, and we report on a search for additional lines in high resolution Chandra 
spectra. No previously unclaimed features are found. We also discuss the significance of lines sets 
reportedly discovered using XMM data for GRB 011211 and GRB 030227. We find that these 
features are likely of modest, though not negligible, significance. 
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I. INTRODUCTION 



■ Multiple luminous X-ray lines have been claimed in spectra taken with XMM of 2 

: GRB afterglows. (GRB 011211, Reeves et al. [9]; GRB 030227, Watson et al. [13]) 

The statistical significance of the GRB 011211 lines has been called into question by 
Rutledge & Sako [11], and we address this question in Section IV. In Butler et al. [1], we 
discuss high resolution spectra from Chandra for the X-ray afterglows to GRB 020813 



O \ and GRB 021004, and we report the discovery of moderately low significance spectral 

\ lines in the case of GRB 020813. Our detection (with an independent instrument) 

■ supports the claimed multiple— a line detections in XMM data, and the high spectral 

resolution facilitates a clearer determination of the line significances (see Section II). 
^ ; There are 3 additional bursts with high resolution Chandra spectra (GRB 991216, 



■ GRB 020405, and GRB 030328), which we have analyzed to search for discrete spectral 

features. Only the data for GRB 991216 yielded a claimed line detection in the literature 
[7]. My collaborators and I analyzed the GRB 030328 data [3, 5, 4] and reported no line 
detections. Null detections for GRB 020405 are reported by Mirabal et al. [6]. 

IL S AND SI LINES FOR GRB 020813? 

In Butler et al. [1], we describe our data reduction and continuum fits for the Chan- 
dra HETGS observation of GRB 020813 (and GRB 021004). We also describe our 
emission line search method which involves successively binning the spectral data by 
factors of 2 in order uncover deviations from the continuum fit. This procedure is sen- 
sitive to resolved emission lines. It allows for quick and easy determination of line 
significances, easily verified by Monte Carlo. One prominent emission line was found 
for GRB 020813. Considering the number of wavelength bins searched, we estimate a 
multiple-trial (i.e. blind search) significance for the line of 3.3c7. If the line is identified 
with the Ka transition in H-like S, a low significance line possibly due to the Ka transi- 



TABLE 1. Lines detected in the spectra of GRB 991216, GRB 020405, and 
GRB 030328 are shown along with the implied redshift (Zx-ray)- This z can be 
compared to that measured in the optical (zopticai)- For GRB 991216, the data do 
not constrain the ionization state of the candidate Fe line. 
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tion in H-like Si can be identified. The significance of the pair turns out to be modestly 
better (3.5a) than the significance of the S line alone [1]. 

Our S line significance determination agrees with the estimate made using the depre- 
cated likelihood ratio test [see, e.g., 8]. To check whether this fact is statistically mean- 
ingful, we apply Monte Carlo integration to establish the true distribution for the log- 
likelihood. We form 10^ simulated data sets using power-law model parameters (with 
Galactic absorption) drawn from the posteriori distribution (the distribution of model 
parameters given the observed data). Each data set is then fitted with this model, then 
with this model plus a Gaussian emission line. The number of iterations yielding larger 
improvements in A;^^ than the observed value (A;(^^ =15.5 for 3 additional degrees of 
freedom) is recorded. To ensure that the parameter space for the emission line is ad- 
equately explored in each Monte Carlo iteration, we use FFTs to determine the most 
significant line-like residual on a fine line centroid wavelength grid (5A = O.OIA) and 
line width grid (dyadic intervals from Ox = 0.02A to = 10.24A). We find that 13 
of 10^ runs yield a larger Ax^ than the observed value, consistent with the significance 
estimates quoted above. We find consistent results independent of whether uniform or 
delta function priors are assumed on the power-law model parameters, indicating that 
the parameters are well constrained by the data. 

III. LINES IN OTHER HIGH RESOLUTION CHANDRA SPECTRA 

We apply the methods which led to the discovery of the lines for GRB 020813 to 
the 4 other high resolution spectra taken with Chandra of GRB X-ray afterglows 
(GRB 991216, GRB 020405, GRB 021004, and GRB 030328). We reduce the data as 
described in Butler et al. [1], with the exception that we jointly fit the gratings data along 
with the Oth-order data, instead of fitting the gratings data alone, in order to best deter- 
mine the continuum fits. The continuum fits and the line searches in the gratings data 
are described in detail in Nat Butler's Ph.D. thesis [2]. The results are shown in Table 1. 
No highly significant, previously unreported lines are detected, and the Fe line reported 
by Piro et al. [7] is robustly detected with a multiple trials significance (conservatively) 
better than 3.7a. This line is apparent in each of the 4 independent Ist-order HETGS 
spectra, and it is also seen in the Oth-order data. It is, in our opinion, the best case for an 
emission line in a GRB X-ray afterglow to date. 

IV. THE SIGNIFICANCE OF THE XMM MULTIPLE -a LINES 

We reduce the EPIC-pn data for GRB 011211 following Reeves et al. [10], finding 537 
net counts in the first 5 ksec of the observation. To search for emission lines, we employ 



the "matched fiker" technique described in Rutledge & Sako [11]. We correct a number 
of minor errors in that work: (1) Rutledge & Sako [11] base their significance estimates 
on model continua determined using fits of sparsely binned (~ 12 counts/bin) data. 
We find that the grouped data bin boundaries are not robust and that modest shifting 
can arise due to minor changes in the source and background selection regions. This 
shifting is sufficient to occasionally wash out the indication of line emission. We choose 
to fit continuum models to the unbinned spectrum. (2) Rutledge & Sako [11] uniformly 
sample from a number of possible continuum models rather than sampling according to 
the posteriori distribution of possible models given the observed data. For example, in 
the case of "x} fitting, our approach would suggest sampling from a model A in frequency 
relative to a model B as exp |— 0.5;f^/;f|}, assuming uniform priors on the model 
parameters. (3) Rutledge & Sako [11] also unjustifiably increase the normalizations on 
their continuum models (relative to the best fit values) in order to force each to yield the 
total number of observed source counts. 
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FIGURE 1. The raw PI spectra for 2 XMM X-ray afterglows are convolved with the line 
response function (i.e. "the matched filter" of Rutledge & Sako [11], dark curve). Significance 
contours are calculated using an MCMC method, which does not rely on binning the data. The 
locations of emission lines claimed by Reeves et al. [9] and Watson et al. [13] are indicated. 

We model the continuum emission as an absorbed power-law, with the absorption 
fixed at the Galactic value as in Reeves et al. [10]. We apply the Markov Chain Monte 
Carlo (MCMC) algorithm described in Appendix B of van Dyke et al. [12] to determine 
the distribution of null model parameters given the data in raw PI bins. The background 
is modelled using a broken power-law with a break at 1.35 keV as in Rutledge & Sako 
[11]. We consider a uniform (i.e. totally non-informative) prior distribution on each 
model parameter for the source and background counts. From 10^ models generated 
in the fashion, we simulate 10^ spectra and apply the matched filter of Rutledge & Sako 
[11] to each. 

The first panel in Figure 1 shows the result of applying this filter to the observed data 
(dark line). The mean result from the simulated spectra is plotted as a dotted line. At the 
locations of the reported S, Si, and Ar emission lines [see, 9], the solid dark curve in 
the first panel of Figure 1 deviates from the mean at > 2o significance. The significance 
contour reached by each feature is plotted, and the individual (single-trial) significances 
of the lines are 3.6a, 2.7a, and 2.4a. Generating an additional 10* model sample 



parameters and reapplying the matched filter to the simulated spectra for each model, 
we count the number of times that a simulated model breaks through each significance 
contour between 0.2 and 2.0 keV. This allows us to determine that the multiple trials 
significances (0)v values in Figure 1) for the emission lines are 2.3<Jn, 1 -Oo^, and OJOn- 
The significance of the line set can then be estimated as 5 times the product of these 
probabilities (i.e. 2Ao), where the factor of 5 takes into account that the data have been 
analyzed in 5 time regions. This is higher than the ~ 1 .0(7 value determined by Rutledge 
& Sako [11], and it is lower than the 3.9(7 initially suggested by Reeves et al. [9]. 

Recently, Watson et al. [13] have claimed the detection of a line set in GRB 030227, 
remarkably similar to the detections claimed for GRB 011211 by Reeves et al. [9]. 
Watson et al. [13] quote a significance determined solely via the likelihood ratio test 
(2.7 —4.4a, depending on the number of degrees of freedom). We reduce the EPIC-pn 
data for GRB 030227 following Watson et al. [13], finding 1593 source counts in the 
final 10 ksec of the observation. In our MCMC analysis, we describe the background 
counts as a power-law. The source counts are modelled as an absorbed power-law as 
in Watson et al. [13]. All model parameter prior distributions are taken as uniform. 
As displayed in the second panel of Figure 1, the claimed Mg, Si, and S lines have 
significances ^ 2.0a (single-trial). This data set is the last of 4 time slices within the 
observation. Performing the multi-trial significance calculation as above, we find the the 
significance of the line set is 1 .3a. Thus, we find that the line emission for GRB 030227 
is less significant than that for GRB 011211. We stress that this significance estimate, 
like that made for GRB 011211 above, is a lower bound on the significance of the line 
emission. If the line centroids can be argued to be constrained by the physics, as argued 
in Watson et al. [13], then the significance would increase. 

V. CONCLUSIONS 

We conservatively estimate the significances for line emission in GRB 991216 (Chan- 
dra), GRB 020813 (Chandra), GRB 011211 (XMM), and GRB 030227 (XMM) as 
> 3.7a, 3.5a, 2.4a, and 1.3a, respectively. We do not find the line emission for 
GRB 011211 (which has galvanized many GRB researchers) to be entirely insignifi- 
cant, as claimed by Rutledge & Sako [11]. However, we find none of the observations to 
be highly significant. Hopefully, early observations with Swift in the coming few years 
will decide conclusively whether this emission is real. 
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